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ABSTRACT 



Observations show that waves are ubiquitous in the solar atmosphere and may play an important role for plasma heating. The study 
of waves in the solar corona is usually based on linear ideal magnetohydrodynamics (MHD) for a fully ionized plasma. However, 
the plasma in the photosphere and the chromosphere is only partially ionized. Here we investigate theoretically the impact of partial 
ionization on MHD wave propagation in cylindrical flux tubes in the two-fluid model. We derive the general dispersion relation that 
takes into account the effects of neutral-ion collisions and the neutral gas pressure. We take the neutral-ion collision frequency as an 
arbitrary parameter. Particular results for transverse kink modes and slow magnetoacoustic modes are shown. We find that the wave 
frequencies only depend on the properties of the ionized fluid when the neutral-ion collision frequency is much lower that the wave 
frequency. For high collision frequencies realistic of the solar atmosphere ions and neutrals behave as a single fluid with an effective 
density corresponding to the sum of densities of both fluids and an effective sound velocity computed as the average of the sound 
velocities of ions and neutrals. The MHD wave frequencies are modified accordingly. The neutral gas pressure can be neglected 
when studying transverse kink waves but it has to be taken into account for a consistent description of slow magnetoacoustic waves. 
The MHD waves are damped due to neutral-ion collisions. The damping is most efficient when the wave frequency and the collision 
frequency are of the same order of magnitude. For high collision frequencies slow magnetoacoustic waves are more efficiently damped 
than transverse kink waves. In addition, we find the presence of cut-offs for certain combinations of parameters that cause the waves 
to become non-propagating. 

Key words. Sun: oscillations — Sun: atmosphere — Sun: magnetic fields — waves — Magnetohydrodynamics (MHD) 
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1. Introduction 

Observations show that magnetohydrodynamic (MHD) waves 
are ubiquitous in the solar at mosphere. Waves are routinely ob- 
served in the corona (e.g.., Mcintosh et al j|201 lb Tomczvk & 
Mcintosh 2009I: iTomczvk et al.ll2007l) , the c hromosp here (e.g.. , 
De Pontieu et al]l2007r.lHe et al]l2009l:lKukhianidze et alJ [2006: 
Oka moto & De Pontieull201 lUZaaarashvili et al.ll2007|). and th e 
photosphere (e.g.,[Fuiimura & Tsuneta H2009l : iJess et al.l 12009). 
In addition, waves and oscillati ons in thin threads of solar promi- 
nences have been reported (e.g., iLin et al.l2007l l2009: Nin g et al.l 
l2009HOkamoto et alj2007h . These recent observations have mo- 
tivated a number of theoretical works aiming to understand the 
properties of the observed waves using the MHD t heory (e.g., 
iPascoe et al JI20 1 2L |20 1 Ol: ISoler et alj|2012al 12011 alibi Terra das 
et al. 120101: 1 Van Doorsselaere et alJl2008t IVerth et alll2010l) . It 
is believed that the observed waves may play an important role 
for the heating of the atmospheric plasma (see, e.g., Cargill & 
de Moortell201 lblDe Pontieu et alJl2007l : lErdelvi & Fedunl2007L 
lMcIntoshetal.1120111) . 

The study of waves in the solar corona is usually based on the 
linear ideal MHD theory for a fullv ionized plasma (see Goed- 
bloed & Poedts l2004l : lPriestlll984 . In this*context. the paper by 



Ed win & Robertsl {l983) on wave propagation in a magnetic flux 
tube has been used as a basic refere nce for subsequent works in 
this field (see also the paper s by, e.g JCally|1986l : lGoossens et al.l 
I2012L 120091 : IWentze]|ll979h . The assumption of full ionization 
is adequate for the coronal plasma. However, the lower tem- 
peratures of the chromosphere and the photosphere cause this 
assumption to become invalid. Therefore, it is necessary to de- 
termine t he impact of partia l ioniz ation on the wave modes dis- 
cussed in Edwin & Roberts ( 1983) in order to perform a realis- 
tic theoretical modeling of wave propagation in partially ionized 
magnetic wave guides. 



There are many papers which have studied the effect of 
neutral-ion collision on MHD waves (see, e.g 



Braginskii 1965. 

De Pontieu et al.ll2001l:lForteza et aljl2007t iKhodachenko et all 



2004120061: Leak e et alJl2005l) . Since the present paper deals 
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with wave propagation in a magnetic cylinder, here we discuss 
those works that studied waves in structured media. In solar 
plasmas the typical frequency of the observed waves is lower 
than the expected value of the neutral-ion collision frequency. 
For this reason the single-fluid theory is usually adopted as a 
first approximation. The single-fluid approximation assumes a 
strong coupling between ions and neutrals, so that bo th fluids 
behave in practice as a single fluid. ISoler et ail (l2009al) studied 
wave propagation in a partially ionized flux tube in the single- 
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fluid approximation. The results of Soler et al. j[200 9a) show 
that the waves discussed in lEdwin & Roberts ( 19831) aredamped 
due to neutral-ion collisions. Alth ough partial ionization may b e 
relevant for plasma heating (see Khomenko & Collados 2012), 
the damping by neutral-ion collisions is in general t oo weak to 
explai n the observed rapid attenuation. In addition, ISoler et al] 
( 2009a) found critical values for the lon gitudinal wavelength tha t 
constrain wave propagation. However. IZaqarashvili et al](l2012l) 
have pointed out that these critical wavelengths may not be phys- 
ical and may be an artifact of the single-fluid approximation. 
This is so because the single-fluid approximation misses impor- 
tant effects when small length scales and/or high frequencies are 
involved. In such cases, the gene ral multifluid description is a 
more suitable approach (see, e.g., IZaqarashvili et al.1l2011l) . In 
a multifluid description no restriction is imposed on the relative 
values of the wave frequency and the neutral-ion collision fre- 
quency. In this more general treatment the mathematical com- 
plexity of the equations is substantially increased compared to 
the single-fluid case. 

A particular form of the multifluid theory is the two-fluid 
theory in which ions and electrons are considered together as an 
ion-electron fluid, i.e., the plasma, while neutrals form another 
fluid that interacts with the plasma by means of collis ions. This 
approach was followed by iKumar & Robertsl d2003l) to study 
surface wav e propagation in a Cartesi an magnetic interface. In 
the study by IKumar & Rober ts (2003) the gas pressure of neu- 
trals was neglected, meaning that the dynamics of neutrals was 
only governed bv the friction force with ions. Recently Soler 
et al. d2012al) used th e two-fluid theory to study resonant Alfven 
waves in flux tubes. I Soler et al] d2012al) focused on transverse 
kink waves, and so they also neglected gas pressure. 

In the present paper we go beyond these previous studies. 
We derive the general dispersion relation for waves in cylindri- 
cal flux tubes in the two-fluid theory. To do so, we consider a 
consistent description of the neutral fluid dynamics that includes 
the effect of neutral gas pressure. This is necessary for a realis- 
tic description of magnetoacoustic waves. The use of cylindri- 
cal geometry and the consideration of neutral pressure are two 
significant improvements with respect to the Cartesian model of 
IKumar & Roberts (2003). The dispersion relation derived here 
is the tw o-fluid generalization of the well-known dispersion re- 
lation of lEdwin & Robertsl (fi 98 3) for a fully ionized single-fluid 
plasma. The collision frequency between ions and neutrals and 
the ionization degree are two important parameters of the model. 
The impact of these two parameters on the waves discussed by 
Edwin & Roberts ( 1983) in the fully ionized case is investigated. 

This paper is organized as follows. Section |2]contains a de- 
scription of the equilibrium configuration and the basic equa- 
tions. In Section[3]we follow a normal mode analysis and derive 
the general dispersion relation for the wave modes. The case in 
which neutral gas pressure is neglected is explored in Section [4] 
both analytically and numeri cally. We compare o ur results with 
the previous findings of Kumar & Roberts (2003). Then, we in- 
corporate the effect of neutral gas pressure in Section Finally, 
we discuss the implications of our results and give our main con- 
clusions in Section[6] 

2. Model and two-fluid equations 

We study waves in a partially ionized medium composed of ions, 
electrons, and neutrals. We use the two-fluid theory in which 
ions and electrons are considered together as an ion-electron 
fluid, while neutrals form another fluid that interact s with the 
ion-electron fluid by means of collisions (see, e.g., iDfaz et a"D 



120121 : ISoler et~ai]l2012allbT: IZaqarashvili et al.ll2011l) . In all the 
following expressions, the subscripts 'ie' and 'n' refer to ion- 
electrons and neutrals, respectively. 

We assume a cylindrically symmetric equilibrium and, for 
convenience, we use cylindrical coordinates, namely r, ip, and 
Z for the radial, azimuthal, and longitudinal coordinates. The 
equilibrium is made of a straight magnetic flux tube of radius R 
embedded in an unbounded environment. The equilibrium mag- 
netic field is straight and constant along the z-direction, namely 
B = B z- Gravity is ignored. As a consequence of the force bal- 
ance condition, the gradients of the equilibrium gas pressures of 
the different species are zero. We also assume a static equilib- 
rium so that there are no equilibrium flows. We restrict ourselves 
to the study of linear perturbations superimposed on the equilib- 
rium state. Hence, the governing equations are linearized (see 
the general expressions in, e.g.. IZaqarashvili et al.ll201ll) . We 
adopt some additional simplifications that enable us to tackle the 
problem of wave propagation analytically. We neglect collisions 
of electrons with neutrals because of the low momentum of elec- 
trons. We also drop from the equations the nonadiabatic mech- 
anisms and the magnetic diffusion terms since our purpose here 
is to determine the impact of neutral-ion collisions only. Thus, 
the set of coupled differential equations governing linear pertur- 
bations from the equilibrium state are 



Pn 



P '~d7 

<9V„ 
dt 

<9b 

dt 

d/7ie 

at 

dt 



-Vpie + - (V X b) X B - p n V ni (Vi - V n ) , (1) 
J" 

-Vp B - p n v ni (V n - Vi) , (2) 

Vx(viXB), (3) 

- r P ie V • Vi, (4) 

-jPrN ■ V„, (5) 



where Vj and v n are the velocities of ions and neutrals, respec- 
tively, pie and p n are the pressure perturbations of ions-electrons 
and neutrals, respectively, b is the magnetic field perturbation, 
Pi and p n are the equilibrium densities of ions and neutrals, re- 
spectively, Pi e and P n are the equilibrium gas pressures of ions- 
electrons and neutrals, respectively, p. is the magnetic permit- 
tivity, y is the adiabatic index, and v n ; is the neutral-ion colli- 
sion frequency. For our following analysis it is useful to define 
the ionization degree as^- = p n /pi- This parameter ranges from 
X — for a fully ionized plasma to^- — > oo for a neutral gas. Also 
for convenience, in the following expressions we replace the ve- 
locities of ions and neutrals by their corresponding Lagrangian 
displacements, £ and £ n , given by 



dt ' 



dt ' 



(6) 



3. Normal modes 

We assume that the equilibrium densities p; and p n are functions 
of r alone, so that the equilibrium is uniform in both azimuthal 
and longitudinal directions. The temperatures of ion-electrons 
and neutrals vary with r accordingly to keep constant the equi- 
librium gas pressures of both fluids. Hence we can write the 
perturbed quantities proportional to exp {wup + ik z z), where m 
and k, are the azimuthal and longitudinal wavenumbers, respec- 
tively. We express the temporal dependence of the perturbations 
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as exp(-/«f), with to the frequency. Now we combine Equa- 
tions (H}-© and arrive at four coupled equations for the radial 
components of the Lagrangian displacements of ions and neu- 
trals, £ r> j e and £ rn respectively, the ion-electron total pressure 
Eulerian perturbation, P' = p K + Bb z /p, and the neutral pres- 
sure Eulerian perturbation, p n , namely 



dP' 
dr 

dp n 
dr 

r dr 

p n c a D 

r or 



'PnVniW^n, 

-ipnVmCjgrM + Pn« (tO + lVm)&,n, 

Vn 



= CiP'+i- 



to + «v, 

,,2 



-C 2 Pn, 



i —PnC n C 2 P' + C 3 Pn, 

tO + lV a i 

with the coefficients T>, C\, C2, and C3 defined as 
2) = Pl (c 2 A +c 2 e )(to 2 -oo 2 A )(tb 2 -to 2 ), 

- (to 2 -to 2 A )(tb 2 -k 2 c 2 e ), 

2 

C 2 = %cl(tf~<^ + ^(c\ + c$(tf~<4)^ 
( 2 

C 3 = 



C 



[(to + zVni) 

D 

to (to + iv n i) 



■Pncl [C 2 + U 2 A (w 2 - W 2 )] 



tO (tO + i'y n i) - C n — t- + k z 



(7) 
(8) 

(9) 

(10) 

(ID 

(12) 
(13) 

(14) 



where u> 2 , a> 2 A , and w 2 are the squares of the modified frequency, 
the Alfven frequency, and the cusp frequency given by 



o?\ \ + 



= kic 



2 ,,2 



kz 



2 2 
2 C a4 



Zc l + C t 



(15) 
(16) 
(17) 



3. 1 . Piece-wise constant equilibrium 

In the present paper we are not interested in studying the res- 
on ant behavio r of the waves. The interested reader is referred 
to ISoler et al.l (l2012aL 12009b) for studies of resonant waves in 
partially ionized plasmas. Instead, here we focus on the modifi- 
cation of the wave frequencies due to neutral-ion collisions. For 
this reason, we avoid the presence of the resonances by choosing 
both p; and Pn to be piece-wise constant functions, namely 



pi,ex 



if r < R, 
if r > R, 



Pn(r) = 



Pn,o 

Pn.ex 



if r < R, 
if r > R, 

(19) 



where R is the radius of the cylinder with constants densities p; o 
and p n embedded in an external environment with constants 
densities p\ , ex and p n ex . This is the density profile adopted by 
lEdwin & Roberts! Fl 98 3) in the fully ionized case. Hereafter the 
indices '0' and 'ex' represent the internal and external media, 
respectively. Note that since the magnetic field and equilibrium 
pressure of both fluids are constant, the rest of equilibrium quan- 
tities are piece-wise constants as well. For simplicity we also 
take v n ; as a constant free parameter. 

For a piece-wise constant equilibrium we can combine Equa- 
tions © and (1 1 Ob to eliminate £ r i e and £ r> „ and obtain two cou- 
pled equations involving P' and p n only, namely 



d 2 P' IdP' I , m- 
or 1 r or 



P' 



d 2 p n 1 d£n_ (jj m 2 . 

dr 2 + r dr r 2 \P* 



<l 2 lPn, 



qlP', 



(20) 
(21) 



with 



[co (to + ixvm) - oj a ) 

(4 + c f e )(^-^) 



cb 2 - 



■k 2 c 2 



~2 2 



K co((o + ixv ni )-co 2 A , 



(22) 



to (to + iVni) - k\c\ 



In addition, c 2 A is the square of the Alfven velocity and c 2 e and 
c\ are the squares of the sound velocities of ion-electrons and 
neutrals, respectively, computed as 



fi 2 



Pn 



(18) 



Equations (|7||-([T0|> are valid for all density profiles in the 
radial direction. For c K = c n — and p„ - they revert to 
the equations discussed in ISoler et al.l d2012ah . Equations (01- 
( llOl l are singular when T> — 0. The positions of the singularities 
are mobile and depend on the frequency, to. The whole set of 
frequencies satisfying T> — at some r form two continua of 
frequ encies called the A lfven and cusp (or slow) continua (see, 
e.g.,[Xppert et aTl l 1974b . A wave whose frequency is within any 
of the two continua is damped by resonant absorption (see, e.g., 
iGoossens et al.ll201ll) . Resonantly dampe d waves in partially 
ionized magnetic cylinders were studied by Sole r et all (120 09b) 
in the single-fluid approximation and by ISoler et al.l (l2012ah in 
the two-fluid formalism. 



01 = 



to + 1 



(23) 



tO + iVnj 



to (to + ftfc) - k 2 c\ k 2 c 2 e {tb 2 - to 2 A ) 



to 2 A to 



to + l\ 



Vni( C 2 +C 2)(^_ w 2) 



XVnitOtO 2 



(c 2 A + c 2 e )(tb 2 -tol) 

(24) 



(c 2 A + cl)(io 2 -oo 2 ) 



(25) 



Equations ( 120b and (121b are the basic equations of this investi- 
gation. They are two coupled Bessel-type differential equations 
representing the coupled behavior of ion-electrons and neutrals 
due to neutral-ion collisions. Note that the terms with q\ and 
q2 are the ones that couple the equations. These terms vanish 
in both fully ionized and fully neutral cases. We need to solve 
Equations (1201 and (|2"TT > to find the general dispersion relation. 
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For later use we also give the expressions of ; e and n in 
terms of P' and p n , namely 



Pi (w 2 -CJ 2 A )\ dr u + iv m dr J ' 



(26) 



■ I 



p n u (oj + iv ni ) (oj + iv ni ) 2 pi (cj2 - w 2 \ ^ 

Vni 1 BP' 

u + ivm pi (w 2 - w 2 ) <3r ' 



dp n 

dr 



(27) 



Note that for v n i = ions and neutrals are decoupled and the 
magnetic field has no influence on the dynamics of neutrals. 

3.2. Strong thermal coupling and characteristic velocities 

Up to here, no restriction on the values of the characteristic 
Alfven and sound velocities has been made. As a result of the 
pressure balance condition at r = R, the six characteristic veloc- 
ities in the equilibrium are related by 

Pi,0 (|4,0 + 4,0 +Xclo) = Pi.ex (| c A,ex + c ie,ex + * C n,ex) ■ ( 28 ) 

Now for simplicity we assume a strong thermal coupling be- 
tween the ionized and neutral fluids, so that both fluids have lo- 
cally the same temperature. However, the temperature can still 
be different in the internal and external media. Using the ideal 
gas law for both species, this requirement has the consequence 
that the local sound velocity of the ionized and neutral fluids are 
related by c 2 = 2c 2 . Hence we can define an effective sound 

J ie n 

velocity of the whole plasma, c s e ff, as 



-s.eff 



2 + X , 



1 + X 2(1 + X ) 16 



(29) 



so that the parameter space of characteristic velocities can be 
reduced from six to four different velocities, namely ca,o, Cs.eff.o, 

CAex, and C s , e ff ; ex- 



3.3. Coupled solutions 

We continue the mathematical study of the normal modes. We 
look for solutions to the coupled Equations ( 120b and (12 1 b . In the 
internal medium, i.e., r < R, physical solutions imply that both 
P' and p n are regular at r — 0. The general solutions of P' and 
p n satisfying this condition in the internal medium are 



P'o = A] J m (k ua r) + A 2 J m (k 2 ,or) , 
p n ,o = _ ^ip ■j^-Jmih.or) 

*1,0 ~~ n,0 
K 2,0 K n,0 



(30) 



(31) 



where J m is the Bessel function of the first kind of order m, A\ 
and A2 are constants, and k\$ and £2,0 are me two possible values 
of the radial wavenumber, k, in the internal medium given by the 
solution of the following equation 



(£ 2 -4,o)(£ 2 -C)-4?,o<72,o = 0. 
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(32) 



Equivalently, in the external medium, i.e., r > R, the solutions 
of P' and p n are of the form 



P' ex = A 3 H™ (k Uex r) + A A H% (*2,„r) , 



Pn,t 



= -A- 



k 2 - k 2 
l,ex n . ex 



tlx 



k l - k 

"-2,ex "-n-ex 



: H^(k 2 , ex r). 



(33) 



(34) 



where H^/ is the Hankel function of the first kind of order m, 
A3 and A4 are constants, and fci^x and k 2fix are the two possi- 
ble values of the radial wavenumber, k, in the external medium 
also given by the solution of Equation d32b but now using ex- 
ternal values for the parameters. The expressions of P' and p n 
for r > R are written generally in terms of Hankel functions 
instead of the usual modified Bessel functions, K m , used to rep- 
resent trapped waves. Thus, our formalism takes into account 
the possibility of wave leakage, i.e., wave radiation in the ex- 
ternal medium. Although we do not explore leaky waves in the 
present work, the obtained dispersion relation will also be valid 
for leaky waves. For ideal trapped waves, the external radial 
wavenumber is purely imaginary and the function ff„ consis- 
tently reverts t o the modified Bessel fu nction of the second kind, 
K m , used by [Ed win & Roberts (1983|). In the partially ionized 
case the radial wavenumber is a complex quantity and the use of 
the function causes us to be very careful when choosing the 
branch of the external radial wavenumber. To avoid nonphysi- 
cal energy propagation coming from infinity we need to select 
the appropriate branch of the external radial wavenumber s o that 
the condition for outgoin g waves is satisfied (see details in lCallvl 
ll986HStenuit et al.lll999l) . 

To understand why there are two different values of the ra- 
dial wavenumber, k, it is instructive to assume a weak coupling 
between the species, i.e., v n i <K \u>\, so that the quadratic terms 
in v„i can be neglected in Equation d32b . The neglected terms 
involve the coupling coefficients q\ and q 2 . Then, the two inde- 
pendent wavenumbers in Equation (132b simplify to 



k\ ~ k[e, k 2 ~ £ n , 



(35) 



where we have dropped the indices '0' or 'ex' because the same 
result is valid in both internal and external media. The two values 
of k for weak coupling reduce to the wavenumbers of the ionized 
and neutral fluids, respectively. For high collision frequencies, 
these two wavenumbers contain terms with q\ and q 2 that couple 
the two fluids. 



3.4. Dispersion relation 

To find the dispersion relation we match the internal solutions to 
the external solutions by means of appropriate boundary condi- 
tions at r — R. The boundary conditions at the interfa ce between 
two m edia in the two-fluid formalism are discussed in Diaz et al. 
(120121) . In the absence of gravity, these boundary conditions re- 
duce to the continuity of P' , p n , £,. j e , and £ r n at r = R. After 
applying the boundary conditions at r — R we find a system of 
four algebraic equations for the constants A\, A 2 , A3, and A4. 
The condition that the system has a non-trivial solution provides 
us with the dispersion relation. The dispersion relation is 



£>,„ (w, k z ) = 0, 



(36) 



with the full expression of T> m given in the AppendixlAl This dis- 
persion relation is the two-fluid generalization of the dispersion 
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relation of Edwin & Roberts! (Il983h for a fully ionized single- 
fluid cylindrical flux tube. In the following sections we discuss 
the corrections due to neutral-ion coll isions on the wave modes 
described by Edwin & Roberti (Il983h in the fully ionized case. 

4. Pressureless neutral fluid 

Before exploring the solutions of the general dispersion relation 
let us consider first the case that neutral pressure is neglected and 
only the gas pressure of the ionized fluid is taken into account. In 
the absence of neutral pressure the only force acting on the neu- 
tral fluid is the friction f orce due to neutral-ion c ollisions. This 
is the situation studied bv lKumar & Roberts! d2003l) in planar ge- 
ometry. We set P a = so that c n = 0. After some algebraic 
manipulations the general dispersion relation simplifies to 



J'JhR) 



KJJc ex R) 



(37) 



with ko and k ex given by 



k 2 - 



(^-<o)(^-^<o) 
Ro + <o)(^ 2 -<o) 

(^-< ex )(^-^U 

( C i,ex+<ex)(^ 2 - W c,ex) 



(38) 
(39) 



In Equation ( 1371 i we have used the modified Bessel function of 
the second kind, K m , instead of the Hankel function of the first 
kind, H^', in order to w rite Equation (|37t in the s ame form as 
the dispersion relation of lEdwin & Robertsl (119831) . Indeed, the 
comparison between Equ ation ( 137b and the dispersion relation of 
Edwin & Roberts ( 1983) reveals that both equations are exactly 
the same if we replace oj 2 by a> 2 in Equation ( 137b . We can take 
advantage of this result to easily study the modifications due to 
neutral-ion collisions of the solutions of the fully ionized case. 



4.1. Analytic approximate study 
4.1.1. Standing waves 

We focus on standing waves so that we fix the longitudinal 
wavenumber, k z , to a real value. Due to temporal damping by 
neutral-ion collisions oj is complex, i.e., oi — wr + icti, where 
<wr and oil are the real and imaginary parts of the frequency, re- 
spectively. Therefore the wave amplitude is damped in time due 
to the exponential factor exp (— \ti>i\f). 

We assume tha t u) — u>q is a trapped solution of the dis- 
persion relation of Edwin & Roberts (1983) in the fully ionized 
case. Here we only co nsider trappe d waves so that uq is real. For 
leaky waves see, e.g., Callv ( 1986). Then Cj 2 — a>^ is a solution 
of Equation (|37| |. Using the expression of u 2 (Equation ( 115b ) we 
expand the Equation u> 2 = u>q as 



+ 1(1 +#)Vn 



0. 



(40) 



The study of the modifications due to neutral-ion collisions of 
the modes of Edwin & Roberts (1983) reduces to the study of 
the solutions of Equation ( 140b . This study is general and is in- 
dependent of the particular mode considered since we have not 
specified what mode corresponds to. All the modes are af- 
fected in the same way by neutral-ion collisions when neutral 
pressure is absent. 



Equation (l40l is a cubic equation, hence it has three solu- 
tions. To determine the nature of the solutions we perform the 
change of variable u> = -is, so that Equation d40b becomes a 
cubic equation in s with all its coefficients real. We compute the 
discriminant, A, of the resulting equation as 

A = -o> 2 [4 (1 + x? < - (x 2 + 20 X - 8) + 4a,*] , (41) 

The discriminant, A, is defined so that (i) Equation ( 140b has one 
purely imaginary zero and two complex zeros when A < 0, (ii) 
Equation ( 140b has a multiple zero and all its zeros are purely 
imaginary when A = 0, and (iii) Equation ( 140b has three distinct 
purely imaginary zeros when A > 0. This criterion points out 
that oscillatory solutions of Equation ( 140b are only possible when 
A < 0. 

To determine when oscillatory solutions are not possible we 
set A = and find the corresponding relation between o>q and v n j 
in terms of x as 
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X 2 +20x-8 X 1/2 (X-V V2 
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(42) 



Since v n i/o>o must be real, Equation (l42b imposes a condition 
on the minimum value of x which allows A = 0. This min- 
imum value is x — 8 and the corresponding critical v n i/o>o is 
Wwo = 1/3 V3 * 0.19. When X < 8 oscillatory solutions are 
always possible. When x > 8 we have to take into account the 
+ and the - signs in Equation d42b . so that Equation ( 142b de- 
fines a range of values of v n i/wo in which oscillatory solutions 
are not pos sible. We call this i nterva l the cut-off region. To our 
knowledge iKulsrud & Pearcd d!969h were the first to report on 
the existence of a cut-of f region of MHD waves in a partially 
ionized plasma, although lKulsrud & Pearcd d 19691) only investi- 
gated Alfven waves in a homogeneous medium. The existence 
of cut-off regions for MH D waves in structured me dia is also 
deduced from the plots of Kumar & Robertsl (120031) but these 
authors did not explore this phenomenon in detail. 

We look for approximate analytic expressions of the fre- 
quency of the oscillatory solutions. Thus we assume that v n i/wo 
is outside the cut-off interval defined by Equation (l42b . meaning 
that Equation ( 140b has two complex (oscillatory) solutions and 
one purely imaginary (evanescent) solution. First we focus on 
the oscillatory solutions. We write u — wr + i<±)\ and insert this 
expression in Equation ( 140b . We assume weak damping and take 
\o>\\ <k |wr|. Hence we neglect terms with o> 2 and higher pow- 
ers. It is crucial for the validity of this approximation that v n j/wo 
is not within or close to the cut-off region in which wr = 0. 
After some algebraic manipulations we derive approximate ex- 
pressions for wr and cjy. For simplicity we omit the intermediate 
steps and give the final expressions, namely 
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X^ + d+xfyl' 

XVni 2 
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(43) 
(44) 



A mode with the same <jj\ but with wr of opposite sign is also a 
solution. When v n ; = 0, so that neutral-ion collisions are absent, 
we find <i>r = coq and a>i = 0. Hence we r ecover the ideal un- 
damped modes of lEdwin & Robertsl (Il983h . Now these modes 
are damped due to neutral-ion collisions and the frequency has 
an imaginary part. In addition, the real part of the frequency is 
also modified. 



Article number, page 5 of 1 1 1 1 



It is instructive to study the limit values of y n j. In the case 



of low collision frequency, 
become 



0»R 

6l>! 



2 



<K co 2 . Equations (l43l and (l44l 

(45) 
(46) 



Thus, the real part of the frequency when <s co^ is just the fre- 
quency found in the fully ionized case (Edwin & Robertsl fT983h 
and it does not depend on the amount of neutrals. The imaginary 
part of the frequency is independent of u>o, meaning that it is the 
same for all the wave modes. On the other hand, for high colli- 
and Equations ( 143 b and ( 144b simplify 



sion frequency, 
to 



0»R 



Oil 



x<4 



(47) 
(48) 



Now the expression of wr involves x m the denominator, so that 
the larger the amount of neutrals, the lower wr compared to the 
frequency in the fully ionized case. The frequency is reduced by 
the factor 1 / Jl + x com pared to the fully ionized value. This is 
the same factor found by Kumar & Roberts (2003) and is equiva- 
lent to replace the ion density by the sum of densities of ions and 
neutrals or, equivalently, to replace p\ by (1 + x)p\- The imagi- 
nary part of the frequency depends now on u>o, meaning that the 
various modes have different damping rates. 

In addition to the oscillatory solutions, Equation (l40l has one 
purely imaginary solution whose approximation is given by 



-2V„ 



»l + (i+x)y 2 ni 



(49) 



The perturbations related to these modes are evanescent in time. 
There exist a different evanescent solution related to each mode, 
coo, of the full y ionized case. Purely i maginary solutions were 
also found by Zaaarashv ili et al.l (1201 ll) when studying waves in 
a partially ionized homogeneous medium. These modes are not 
present in the fully ionized case. These evanescent perturbations 
may be relevant during the excitation of the waves, since part 
of the energy used to excite the waves may go to these modes 
instead of being used to excite the oscillatory solutions. This 
could be investigated by going beyond the present normal mode 
analysis and solving the initial value problem. 

4.1.2. Propagating waves 

Now we turn to propagating waves. In the fully ionized case 
the study of propagating waves is equivalent to that of standing 
waves. Here we shall see that neutral-ion collisions break this 
equivalence and propagating waves are worth being studied sep- 
arately. Propagating waves were not investigated bv Kumar & 
Roberts ( 2003). For propagating waves, we fix the frequency, u>, 
to a real value and solve the dispersion relation for the complex 
k z , i.e., k z - k Z R + ik z j, where £r and k\ are the real and imag- 
inary parts of k z , respectively. Therefore the wave amplitude is 
damped in space due to the exponential factor exp (— k z jz). 

Equation ( I37b also holds for propagating waves. We assume 
that k z - fc 7 .n, with k 7 .n real, is a so lution of the ideal dispersion 
relation of Edwin & Roberts! (l!983l) . Then, 
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is a solution of Equation ( I37b following the same argument as 
explained in the case of standing waves. We write k z = k z ^ + 
ik z j and insert this expression in Equation (15 0) . We obtain exact 
expressions for k z $_ and k z j, namely 
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(51) 



(52) 



The same values of k z ^ and k z j but with opposite signs are also 
solutions. In principle, Equation ( 15 lb allows two different values 
of Lr because of the + sign. However the solution with the - 
sign is not physical since it corresponds to £ z r imaginary, which 
is an obvious contradiction. For this reason we discard the solu- 
tion with the - sign and take the + sign in Equation (IBTI ). Now 
we realize that 



2 2 2 

X v tor 



« 1. 



(co 2 + (l + x)v 2 ^ 
Thus we can approximate r and k z j as 



(53) 
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(55) 



These expressions are equivalent to Equations (143) and (1441 ob- 
tained for standing waves. However, contrary to Equations (I43l > 
and (l44l . Equations (l54l and (155b are valid for all values of the 
ratio v n j/w. There is no cut-off region for propagating waves 
because L r never vanishes. There are no purely imaginary so- 
lutions of k z , i.e., all the solutions have always an oscillatory 
behavior in z. This is an important difference with the case of 
standing waves. 

4.2. Comparison with numerical results 

Here we numerically solve the dispersion relation (Equa- 
tion (137b). For simplicity we focus on standing waves and show 
the results for transverse kink (m = 1) waves only, although we 
have checked that equivalent results are found in the case of other 
modes as, e.g., slow magnetoacoustic modes. In the fully ionized 
case, the frequency of the transverse kink waves in the thin tube 
(TT) limit, i.e., k z R <K 1, is o> — oj^, with co^ the kink frequency 
given by 



u k = k. 



/ Pi.0ci,0 + fe* c lex 
I Pi,0 + Pi.ex 



(56) 



Figure Q] displays the real and imaginary parts of the fre- 
quency of the kink mode as functions of the ratio v n ;/wk for a 
particular choice of parameters given in the caption of the figure. 
We compare the numerical results with the analytic approxima- 
tions. In Figure Q] we have considered^- < 8 so that no cut-off 
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R. Soler, A. J. Diaz, J. L. Ballester, & M. Goossens: Waves in a two-fluid flux tube 




Fig. 1. (a) w R /w k and (b) c^/wjc versus v n i/<jJ k for the transverse 
kink (m = 1) mode with k : R = 0.1, Pi,o/Pi,ex = 3, Ci e ,o/cA,o = 0.2, and 
X = 4. The solid line is the result obtained by numerically solving the 
dispersion relation in the absence of neutral pressure (Equation (137b). 
The symbols are the approximate analytic results in the weak damping 
approximation (Equations d43t and <44t ). 



region is present. Regarding the real part of the frequency (Fig- 
ureQJ), we obtain that for v n ;/aJk ^ 1 the results are independent 
of the ionization degree, with wr * Wk as in the fully ionized 
case. When the ratio v ni /a)^ increases, plasma and neutrals are 
more coupled and the ionization degree becomes a relevant pa- 
rameter, so that decreases until de value wr * co^/ -y/1 + x 
is reached. This behavior is consistent with the analytic Equa- 
tion (03]) although the details of the transition are not fully cap- 
tured by the approximation. Regarding the imaginary part of the 
frequency (Figure QJ)), we find that a>i tends to zero in the lim- 
its Vni/wk <k 1 and v n ;/wk s> 1. The damping is most efficient 
when v n ; and «k are of the same order of magnitude, approx- 
imately. This result is also consistent with the analytic Equa- 
tion d44l) . although the approximation underestimates the actual 
damping rate when co\ is minimal. This is so because Equa- 
tion ( 144b was derived in the weak damping approximation, i.e., 
|a»i[ <s |o»r|, while \a>\\ and |wr| are of the same order when a>\ is 
minimal, meaning that the damping is strong. 

Now we take \ > 8 and repeat the previous computations. 
These results are shown in Figure |2] We notice the presence 
of a cut-off region for a certain range of v n i/wk. In this cut- 
off region wr = 0. The location of the cut-off region agrees 
very well with the range given by Equation (142) . This zone 
exists for relatively small v n j/wk- In the solar atmosphere the 
expected value of the collision frequency (see, e.g., De Pontieu 
et al. l2001l , Figures 1 and 2 ) is much higher than th e frequency 
of the observed waves (e.g. JPe Pontieu et al. 2007; Okamoto & 



Fig. 2. Same as Figure[T]but with^ = 20. 
the cut-off r egion given bv E quation (142b . 



The shaded area denotes 



De Pontieu 120 lC The minimum value of the neutral-ion colli- 
sion frequency in the solar atmosphere is of the order of 10 Hz, 
while the dominant f requency in the observations o f chromo- 
spheric kink waves by |Okamoto & De Pontiiu| (1201 lh is around 
22 mHz. Therefore, the cut-off region found here is not in the 
range of v n i/wk consistent with solar atmospheric parameters and 
observed wave frequencies. However, in other astrophysical sit- 
uations the presence of the cut-off region may be relevant (see 
iKulsrud & Pearcdl 19691) . As expected, the analytic approxima- 
tions (Equations ( 143b and (144b ) completely miss the presence of 
the cut-off region, although they are reasonably good far from 
the location of the cut-off. 

It is worth mentioning again that the effect of resonant ab- 
sorption on the damping of transverse kink waves is absent due 
to the use of a piece-wise constant density profile (Goossens 
et al. 120021: iRuderman & Robertsll2002l) . ISoler et al.l (l2012ah 
showed that damping by resonant absorption is more efficient 
than damping by neutral-ion collisions when v n i/wk » 1. 
Hence, the actual damping of transverse kink waves would be 
stronger than the damping shown he re if resonant abs orption 
were taken into account (see details in lSoler et al.ll2012ab . 



5. Role of neutral pressure 

We incorporate the effect of neutral pressure. We anticipate 
that neutral pressure is relevant for slow magnetoacoustic waves, 
while its effect on transverse kink waves is minor. Goossens 
et al. (120121 12009) showed that transverse kink waves in flux 
tubes have the typical properties of surface Alfven (or Alfvenic) 
waves. In thin tubes, i.e., k z R <K 1, the main restoring force 
of these waves is magnetic tension and the gas pressure force 
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is negligible. We have computed the kink mode frequency in 
the presence of neutral pressure and have found no significant 
differences with the result in the absence of neutral pressure. 
Hence, here we focus on the study of slow magnetoacoustic 
waves and explore the modification of the slow mode frequency 
due to neutral-ion collisions when neu tral pressure is taken int o 
account. In the fully ionized case (Edwi n & Robertsl fl983). 
slow magnetoacoustic waves in the TT limit have frequencies 
a) a o) c ,o, with o) Ci o the internal cusp frequency given by 



oj cfl = k z 



2 2 
C A,0 C ie,0 

C A,0 T °ie,0 



(57) 



This result is independent of the azimuthal wavenumber, m. In 
a low-/? plasma, i.e., when the magnetic pressure is much more 
important than the gas pressure, <y c ,o ~ k z Ci e $- Due to the com- 
plexity of the general dispersion relation (Equation (l36l ) it is not 
possible to obtain simple analytic approximations of the slow 
mode frequency. For this reason we perform this investigation in 
a numerical way. 

Figure [3] shows the numerically obtained real and imagi- 
nary parts of the slow mode frequency as functions of the ra- 
tio v n i/w Cj o ior x < 8. The rest of parameters are given in the 
caption of the figure. For comparison, we overplot the results 
in the absence of neutral pressure which are obtained by solving 
Equation ( 137l i with the same parameters. First we discuss the be- 
havior of the wr (Figure |3^). When y n i/a»c,o « 1 we recover the 
result in the fully ionized case, i.e., o>r ~ w c q. When v n i/ai Ci o 
increases, the real part of the slow mode frequency decreases 
until a plateau is reached at the value wr ss <y Ci o/ -\A + X- This 
is equivalent to perform the replacements c\ — > c 2 A /( 1 + x) an d 
c? — > cf e /(l +x) m Equation (157b and is the same behavior as that 
found in the absence of neutral pressure, so that both results are 
superimposed in Figure [3^. However, as the ratio Vni/a> c> o con- 
tinues to increase, the real part of the frequency rises again until 
a second plateau is finally reached. The presence of this second 
plateau of wr for large collision frequencies owes its existence to 
the effect of neutral pressure and is absent when neutral pressure 
is neglected. 

Let us analyze in detail the second plateau of wr at high colli- 
sion frequencies. Since we have no simple analytic expressions 
for the slow mode frequency we have performed a parameter 
study. We have considered various values of c; e and have numer- 
ically computed the slow mode ojr at the second plateau. The 
result of this parameter study (not shown here for simplicity) 
points out that the value of wr at the second plateau is, approxi- 
mately, 



Wr ~ COcfl 
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This expression can be obtained from Equation (157) by perform- 
ing the replacements c 2 A — > c 2 A /{ 1 + x) and c 2 



i.eff' 



with the 

expression of the effective sound velocity, c Sje ff, given in Equa- 

k z c le o so that 



(59) 



tion ((29). In a low-/? plasma, c A » cf and a> c 
Equation d5"8l simplifies to 



wr » k z c ie 



2+X 
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Hence, for practical purposes we obtain that at high collision 
frequencies the sound velocity of the ionized fluid has to be re- 
placed by the effective sound velocity in the expressions of the 
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Fig. 3. (a) a> R /(i> c0 and (b) LL>i/a> c0 versus v ni /a> c o for the slow mag- 
netoacoustic mode with^- = 4 and the rest of parameters the same as in 
FigureQ] The solid line is the full result taking neutral pressure into ac- 
count, while the dashed line is the result when neutral pressure is absent. 



frequency. The effective sound velocity takes into account the 
sound velocity of neutrals in addition to the sound velocity of 
ions. On the contrary, for intermediate collision frequencies it is 
enough to replace c K by Ci e / yT+x, so that the sound velocity 
of neutrals can be ignored. 

We turn to the imaginary part of the frequency (Figure |3j>). 
As for the real part of the frequency, there are striking di fferences 
between the results with and without neutral pressure. The full 
result shows the presence of two different minima of oj\, whereas 
only the first minimum is present in the absence of neutral pres- 
sure. The second minimum takes place in the range of high col- 
lision frequencies compared to the ideal cusp frequency. As a 
consequence of the presence of this second minimum for high 
collision frequencies, the efficiency of neutral-ion collisions for 
the damping of the slow mode is larger than for transverse kink 
modes when v n j approaches realistic values. The stronger damp- 
ing of the slow modes due to neutral-ion collisions compared to 
the weak damping of the tran sverse kink modes was already ob- 
tained by Sole r et all ([2009a) in the single-fluid approximation. 

Now we consider x > 8 and repeat the previous computa- 
tions (see Figure |4] for x - 20). As in the case without neutral 
pressure, we notice the presence of a cut-off region that agrees 
well with the range given by Equation (1421 . The existence of 
the cut-off region is not affected by neutral pressure, i.e., the lo- 
cation of the cut-off region in the same when neutral pressure 
is neglected. For the set of parameters considered in Figure |4] 
the cut-off region appears in the vicinity of v n i/w Ci o ~ 10 -1 . We 
increase the ionization fraction to x — 100 (Figure [5}. In addi- 
tion to the cut-off region described above, now we find a second 
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Fig. 4. Same as Figure[3]but with^- = 20. The shaded area denotes 
the cut-off region given by Equation d42t . 

cut-off region. The second cut-off region appears for relatively 
high collision frequencies, v n i/a» Cj o ~ 10 1 for the set of param- 
eters used in Figure [5] This new cut-off region is only present 
when neutral pressure is taken into account. The result in the ab- 
sence of neutral pressure completely misses this second cut-off 
region. The existence of a second cut-off region seems a result 
of slow modes only, since it is not found in the case of transverse 
kink modes. Unfortunately, unlike the first cut-off region we do 
not have an analytic expression that gives us the location of the 
second cut-off region. Instead, out numerical study informs us 
that the second cut-off appears when^ > 24, although this value 
may be affected by the particular choice of parameters used in 
the model. 

As for the transverse kink modes, it is worth mentioning 
that some mechanisms not considered in the present analysis 
may produce damping of the slow mode comparable or even 
stronger than that due to neutral-ion collisions. This is the case 
of non-adiabatic mecha nisms, e.g., thermal conduction and ra- 
diativ e losses (see, e.g. JDe Moortel & Hoodll2003HSoler et alJ 
2008). Partial ionization increases the efficiency of thermal con- 
ducti on due to the additional contribution of neutral conduction 
(see lForteza et al Jl2008t ISoler et al.ll2010h . so that the combined 
effect of thermal conduction and partial ionization produces a 
very efficient damping of the slow mode (see, e.g., Khodachenko 
et al. l200dlSoleil2010h . 

6. Discussion 

In this paper we have explored the impact of partial ionization on 
the properties of the MHD waves in a cylindrical magnetic flux 
tube using the two-fluid formalism. Unlike previous works (e.g., 



Goossens: Waves in a two-fluid flux tube 
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Fig. 5. Same as Figure|4]but with^ = 100. 



Kumar & Roberts 2003; [Soier et al.ll2012ah we have considered 
a consistent description of the dynamics of the neutral fluid that 
takes gas pressure into account. We have derived the dispersion 
relation for the wave modes that is the two-flui d gene ralization 
of the dispersion relation of lEdwin & Robertsl (119831) . Rather 
than performing particular applications, we have considered the 
neutral-ion collision frequency as an arbitrary parameter. The 
modifications of the wave frequencies due to neutral-ion colli- 
sions have been explored. 

First, we have investigated the case that neglects neutral 
pressure. This case was previo usly investigated for su rface 
waves in a Cartesian interface by iKumar & R oberts (2003). For 
Vni/M « 1 we have recovered the frequencies in fully ionized 
plasma (lEdwin & Robertslfl983l As the ratio v n j/|w| increases, 
both ion-electrons and neutrals become more and more coupled. 
For Vni/M ^ 1 ion-electrons and neutrals behave as a single 
fluid. Compared to the fully ionized case, when v m 7M » 1 the 
wave frequencies are lower and depend on the ionization ratio, 
X- The frequencies are reduced by the factor (\ + I) -1 ' 2 com- 
pared to their values in the fully ionized case. This result is the 
same for all wave modes and is equivalent to replace the den- 
sity of ions by the sum of densities of ions and neutrals. The 
same conclusion was reached in the previous work bv Kumar & 
Roberts (2003). Concerning the damping by neutral-ion colli- 
sions, we find that damping is weak when Vm/M <s 1 unless the 
plasma is very weakly ionized, i.e., ^ » 1. Damping is most 
efficient when the wave frequency is of the same order of mag- 
nitude as the collision frequency. Damping is again weak when 
Vm/M ^ 1- This last result is always true regardless the value 
of y. Again, our results a gree well with the previous findings of 
IKumar & Robertsl (l2003h . 
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Then we have incorporated the effect of neutral gas pressure. 
Since transverse kink waves are largely insensitive to the sound 
velocity, we have obtained the same results as in the absence of 
gas pressure. However, the consideration of neutral pressure has 
dramatic consequences for slow magnetoacoustic modes. For 
high collision frequencies the slow mode frequency is larger than 
the value obtained in the absence of neutral pressure. This is so 
because the slow mode frequency depends on an effective sound 
velocity that corresponds to the weighted average of the sound 
velocities of ions-electrons and neutrals. When ^- » 1 the ef- 
fective sound velocity tends to the sound velocity of neutrals. 
Concerning the imaginary part of the frequency, the results with 
and without neutral pressure also show significant differences for 
high collision frequencies. The imaginary part of the slow mode 
frequency has an additional minimum at high collision frequen- 
cies when neutral pressure is included. This causes the slow 
mode to be more efficiently damped by collisions when neutral 
pressure is taken into account. 

Here we go back to the discussion given in the Introduc- 
tion about the necessity of using the multi-fluid theory instead 
of the simpler single-fluid MHD approximation. In the solar 
atmosphere the expected value of th e neutral-ion collision fre- 
quency (see iDe Pontieu et al. 2001, Figures 1 and 2) is much 
higher than the frequency of the observed waves (e.g., De Pon- 
tieu et al. 2007; Okamo to & De Pontieull201 lh . In such a case, 
the results of this paper point out that, for practical purposes, it 
is enough to use the single-fluid MHD approximation, but tak- 
ing into account the following simple recipe to adapt the results 
of fully ionized models to the partially ionized case. First of 
all, the ion density has to be replaced by the total density, i.e., 
the sum of densities of ions and neutrals. Second, the effective 
sound velocity should be used instead of the sound velocity of 
the ionized fluid. This effective sound velocity is the weighted 
average of the sound velocities of ions-electrons and neutrals. 
This approach is appropriate for those works which are not in- 
terested in the details of the interaction between ions and neu- 
trals, but only on the result of this interaction on the wave fre- 
quencies. Hence, this recipe may be useful to perform magneto- 
seismology of chromospheric spicules and other partially ion- 
ized structures as prominence threads. For example, in their 
seismological analysis of transverse kink waves in spicules Verth 
et al. d2011l ) considered a definition of the kink velocity that in- 
corporates the density of neutrals. Although our recipe provides 
a good approximation to obtain the wave frequencies, it misses 
the effect of damping due to neutral-ion collisions, which might 
have an impact on the seismological estimates using the ampli- 
tude of the waves. Taking into account the range of observed 
wave frequencies, the damping of kink waves due to neutral-ion 
collisions may be n egligible compared to other effects as, e.g, 
resonant absorption (Soler et al. 2012aj). However, the damping 
of slow magnetoacoustic waves may be important even for high 
collision frequencies. Other damping effects as, e.g., resonant 
absorption and non-adiabatic mechanisms, should be considered 
in addition to neutral-ion collisions to explain the damping of 
MHD waves in the solar atmosphere (see, e.g., Khodachenko 
et al. 120061: ISoler et al.ll2012al) 

On the other hand, a relevant result obtained in this paper is 
the p resence of cut-offs for ce rtain combinations of parameters 
(see iKulsrud & PearceH l969). This result cannot be obtained 
in the single-fluid MHD approximation. At the cut-offs the real 
part of the frequency vanishes, meaning that perturbations are 
evanescent in time instead of oscillatory. For transverse kink 
waves we find a cut-off region for low values of v n ;/wk when 
X > 8. For slow magnetoacoustic waves there are two cut-off re- 



gions. The first one is the same as that found for transverse kink 
waves, whereas the second one is only present in the case of slow 
modes when neutral pressure is considered. The second cut-off 
region takes place when^ > 24 and relatively high collision fre- 
quencies. The relevance of these cut-offs for wave propagation 
in the solar atmosphere should be investigated in forthcoming 
works. 

Finally, we are aware that the model used in this paper is 
simple and misses effects that may be of importance in real- 
ity. The model should be improved in the future by consider- 
ing additional ingredients not included in the present analysis. 
Among these effects, the variation of physical parameters, e.g., 
density, ionization degree, etc., along the magnetic wave guide, 
ionization and recombination processes, and the influence of a 
dynamic background are worth being explored in the near fu- 
ture. 
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Appendix A: Expression of the dispersion relation 
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The dispersion relation is 2)„, (to, k z ) = 0, with D m (to, k z ) given 
by the solution of the following determinant, 
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